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The ability of a recently developed square-integrable discrete basis to represent the properties 
of the continuum of a two-body system is investigated. The basis is obtained performing a simple 
analytic local scale transformation to the harmonic oscillator basis. Scattering phase-shifts and the 
electric transition probabilities B(E1) and B(E2) have been evaluated for several potentials using 
the proposed basis. Both quantities are found to be in excellent agreement with the exact values 
calculated from the true scattering states. The basis has been applied to describe the projectile 
continuum in the 6 He scattering by 12 C and 208 Pb targets at 240 MeV/nucleon and the n Be 
scattering by 12 C at 67 MeV/nucleon. The calculated breakup differential cross sections are found 
to be in very good agreement with the available experimental data for these reactions. 
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I. INTRODUCTION 

Quantum collisions involving weakly bound systems are known to be influenced by the coupling to the unbound 
states. In nuclear collisions, this was first evidenced in the pioneering study by Johnson and Soper who recog- 
nized the relevance of the breakup channels to understand deuteron induced reactions. In their survey, the deuteron 
continuum was represented by a single s-state. Later developments by Rawitscher [2J and Austern [3[ helped to 
introduce a more realistic representation of the continuum, leading to the development of the Continuum-Discretized 
Coupled-Channels (CDCC) method. This method reduces the many-body problem to an effective three-body prob- 
lem and expands the full three-body wave function in a selected set of continuum wave functions of a given pair 
subsystem Hamiltonian. Projection of the Schrodingcr equation onto the selected internal states gives rise to a set 
of coupled differential equations. The practical implementation of the method requires a discretization procedure, 
i.e., an approximation of the two-body continuum spectrum by a finite and discrete representation. Although not 
strictly necessary, it results numerically advantageous to use for this representation a set of L 2 (i.e., square- integrable) 
functions. The standard method of continuum discretization consists on dividing the continuum into a set of energy 
or momentum intervals. For each interval, or bin, a representative wave function is constructed by superposition of 
the scattering states within the interval (the average method) . 

An alternative to the discretization method based on bins is the pseudo-state (PS) method, in which the wave 
functions describing the internal motion of the projectile are obtained as the eigenstates of the projectile Hamiltonian 
in a truncated basis of square-integrable functions. A variety of PS basis have been proposed in the literature for 
two-body continuum discretization and, more recently, also for the three-body continuum 

In a recent work Q, we proposed a PS method based on a Local Scale Transformation (LST) of the Harmonic 
Oscillator (HO) basis. The LST, adopted from a previous work of Karataglidis et al. [ll|, is such that it transforms 
the Gaussian asymptotic behavior into an exponential form, thus ensuring the correct asymptotic behavior for the 
bound wave functions. The accuracy of this THO basis was tested for several reactions induced by deuteron and halo 
nuclei, showing an excellent agreement with the standard binning method, and an improved convergence rate. 

Due to their vanishing asymptotic behavior, it is not obvious that genuine continuum properties, such as the 
scattering phase-shifts, can be well described using square-integrable states. However, in this case one can make use 
of integral formulas, which require only the wave function within a finite region. Several prescriptions have been 
proposed in the literature to extract the phase-shifts from continuum-discretized states [12Hl4j . In this work, we will 
make use of the stabilization method of Hazi and Taylor [l2|, [H, EH to show that the THO basis reproduces very 
well the exact phase-shifts. As an additional test of the quality of the THO basis, we will calculate several transition 
probabilities and their associated sum rules. Finally, we will apply the THO basis to calculate the breakup of the 
reactions 6 He+ 12 C and 6 He+ 208 Pb at 240 MeV/nucleon and n Be+ 12 C at 67 MeV/nucleon, making use of the CDCC 
formalism. The calculations will be compared with existing experimental data for these reactions. 

The work is structured as follows. In Section |TT] we review the THO method based on the parametric LST. In 
Section HTT1 we introduce an integral formula suitable for the calculation of scattering phase-shifts with PS functions. 
In Section IIVI we recall some useful formulae to evaluate the dipole and quadrupole transition probabilities from the 
scattering states and from the pseudo-states. In Section [V] these formulae are applied to study the continuum of the 
deuteron and 6 He nuclei. In Section IVT1 the method is applied to the scattering of 6 He by 12 C and 208 Pb at 240 MeV 
per nucleon and 11 Be+ 12 C at 67 MeV/nucleon. Finally, in Section IVlII we summarize the main results of this work. 



II. THE ANALYTIC LST 

In this section, we briefly review the features of the PS basis used in this work. This basis was originally developed 
in [TTI ] to describe the single-particle orbitals within a mean- field approach. In a later work [f| , we adopted this method 
to discretize the continuum of a two-body system within the context of the Continuum-Discretized Coupled-Channels 
method. 

The starting point is the HO basis in angular momentum representation. The radial part of the n-th HO function 
for a given partial wave I is here denoted 4>nf(s). These functions are orthogonal and constitute a complete set and, 
therefore, they can be used to expand the eigenstates (bound and unbound) of an arbitrary potential. For a finite 
well, the bound state wave functions decay exponentially at large distances and hence the HO basis does not provide 
a suitable representation due to its Gaussian asymptotic form. A possible approach to overcome this limitation, while 
retaining the appealing properties of the HO basis, is to perform a local scale transformation (LST) that converts the 
Gaussian behavior into an exponential one [13, Ell ■ This gives rise to the so called Transformed Harmonic Oscillator 
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(THO) basis. We will denote the radial part of these basis states as: 

<?°W = /f<?[ S (r)]. (1) 

Note that, by construction, the family of functions <f>ni( r ) are orthogonal and constitute a complete set with the 
following normalization: 

poo 

dr|<f°(r)| a = l. (2) 



Moreover, they decay exponentially at large distances, thus ensuring the correct asymptotic behavior for the bound 
wave functions. In practical calculations a finite set of functions ([1]) is retained, and the internal Hamiltonian of 
the projectile is diagonalized in this truncated basis with N states, giving rise to a set of eigenvalues and their 
associated eigenfunctions, denoted respectively {e„} and {vl^O")} ( n — 1, ...,JV). As the basis size is increased, 
those eigenstates with negative energy will tend to the exact bound states of the system, while eigenstates with 
positive eigenvalues can be regarded as a finite representation of the unbound states. 

With the criterion given above, the LST is indeed not unique. In Ref. [l9j the LST was defined in such a way 
that the first HO state is exactly transformed into the exact ground state wave function. Therefore, by construction, 
this wave function is exactly recovered for any arbitrary size of the basis. In a more recent work [6| we adopted the 
parametric form of Karataglidis et al. 



s(r) 



V2b 



(3) 



that depends on the parameters m, 7 and the oscillator length b. Following |6j], the oscillator length b is treated 
as a variational parameter used to minimize the ground state energy. Asymptotically, the function s(r) behaves 
as s(r) ~ j^JT^ and hence the functions obtained by applying this LST to the HO basis behave at large distances 
as exp(— 7 2 r/2& 2 ). Therefore, the ratio 7/6 can be related to an effective linear momentum, fc e ff = 7 2 /26 2 , which 
governs the asymptotic behavior of the THO functions; as the ratio 7/6 increases, the radial extension of the basis 
decreases and, consequently, the eigenvalues obtained upon diagonalization of the Hamiltonian in the THO basis tend 
to concentrate at higher energies. Therefore, 7/6 determines the density of PS as a function of the excitation energy. 
This property was used in [6[ to determine a suitable value for the ratio 7/6 in scattering calculations. For a more 
quantitative measurement of the density of states we define the magnitude: 

JV 

P (N \k) = J2(Mk)\<P ( r3)' ( 4 ) 

71=1 

with \ipi(k)} denoting the scattering wave function for a momentum k. 

With this definition the integral of the density with respect to the momentum is just the number of basis states, 
i.e. 

dkp {N \k) = N : (5) 

regardless of the choice of the parameters of the LST. 

In all the calculations presented in this work, the power m is just taken as m = 4. This was one of the choices done 
in Ref. and, in fact, the authors of that work found a very weak dependence of the results on this parameter. 

III. EXTRACTING THE PHASE-SHIFTS FROM THE THO BASIS 

The properties of the continuum states are completely determined by the phase-shifts. In a two-body problem, 
the phase-shifts are readily obtained from the asymptotics of the radial part of the wave function. Ignoring spins for 
simplicity, the radial part corresponding to a partial wave I can be written at large distances as: 

<pt(k,r) -> ^[cQB8 t (k)F t (kr) + sm8i(k)Gt(kr)] , (6) 
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where F( and Gi are the regular and irregular Coulomb functions. If the potential is real, the functions tpe as well as 
the phase-shifts S# are also real. 

Equation ([6]) can not be applied to PS functions to extract the phase-shifts, because these functions vanish asymp- 
totically. However, the phase-shifts can be also obtained from integral expressions, which require only the interior 
part of the wave functions. Here, we make use of the integral formula proposed by Hazi and Taylor [l2l Il6j . who 
applied this formula to extract the phase-shifts in a one-dimensional scattering problem using a harmonic oscillator 
representation. We have generalized this formula to three-dimensional cases. The formula so obtained reads 

. r /i\ f °°Mk,r)[E-H]f(r)F e (kr)dr 

tandi (k) — — — - — — — - — . (7) 

fA ' / °° Mk,r)[E - H]f(r)G e (kr)dr K) 

This formula can be derived following the same arguments outlined in Ref. [l2[ for the one-dimensional case. We note 
that, if the exact wave functions are used for tpi(k,r), this expression becomes an alternative to Eq. (j6|) to calculate 
the exact phase-shifts. The function f(r) appearing in Eq. ([7]) verifies the following properties: 

/(r)-^l; /(0) = /'(0)=0. (8) 



Following 12j |. we adopt the explicit form f(r) — 1 — exp(— /3r 2 ), with j3 > 0. The aim of this function f(r) is to avoid 
evaluating the function Gt(kr) at the origin, where it becomes singular. Therefore the parameter /3 should be small 
enough to make f(r) f=s for distances of the order of the nuclear range. In the cases studied here, we have chosen 
^ = 0.01fm- 2 . 

IV. ELECTRIC TRANSITION PROBABILITIES IN THE PS BASIS 

The accuracy of the PS basis to represent the continuum can be studied by comparing the ground-state to continuum 
transition probability due to a given operator. Here we consider the important case of the electric dissociation of the 
initial nucleus a into the fragments b + c. This involves a matrix element between a bound state (typically the ground 
state) and the continuum states. 

The electric transition probability between two bound states \(£is)ji) and \(£fs)jf) (assumed here to be normalized 
to unity) is given by the reduced matrix element 

B{EX;i -+ /) = ||±1 W fS ) 3f \\M{E\)\\{l lS ) 3l )\ 2 , (9) 



where M. is the multipole operator: 



M{EXfi) = Z^er A Y x „(f), (10) 



with the effective charge 



\rrib + m c J \ m b + rr 



In the case of a transition to a continuum of states, \k(£fs)jf), the preceding definition is replaced by (see e.g. |20|): 

dB(EX) = 2j> + 1 Mbcfc 
de 2ji + 1 (2ir) s h 2 

x Kfc(^/*)i/| 1^(^)11^*)^) | 2 , (12) 

with k — y/2fl^e/h. Note that the extra factor appearing in Eq. (IT21 with respect to Eq. (0) is consistent with the 
convention (k(£s)j\k'(£s)j) = 5{k — k') and the asymptotic behavior of Eq. ©. 

In the calculations presented in this work, we will ignore for simplicity the internal spins of the clusters and hence 
s = 0, ji = £i, jf — if. In addition, we will consider only transitions from the ground state, so (r\(£is)ji) = y g . s .(r) 
(where the index £j is omitted for shortness). This reduces expression (1121) to: 



(13) 
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The reduced matrix element is given by: 



47T 



(^ f (AOI|M(i?A)||^. s .) = ^Z eS eD^ if Rl ef (k) 



with a geometric factor [20| and R^. if (k) the radial integr 



al 



(14) 



(15) 



Using a finite basis, one may calculate only discrete values for the transition probability. According to Eq. ([9]), the 
B(EX) between the ground state and the n-th PS is given by: 



B<~ N \EX;g.s. -> 



2^/ + l 
2L + 1 



(^l\\M{E\)\\y g ., 



(16) 



In order to relate this discrete representation with the continuous distribution f| 13[) one may use the simple approxi- 
mation: 



dB(EX) 



ds 



-LsW(£A;g.s.^n), 



(17) 



where A„ = (e n +i — £n-i)/2 is an estimate for the energy width of the n-th PS. This expression provides the B(EX) 
values only for the PS eigenvalues s n . 

Alternatively, one may derive a continuous approximation to (|13[) by introducing the identity in the truncated PS 
basis, i.e.: 



N 



(18) 



For TV — > oo this expression tends to the exact identity operator for the Hilbert space spanned by the eigenfunctions of 
the considered Hamiltonian. By inserting (|18j) into the exact expression (|13j) we obtain the approximated continuous 
distribution: 

dB(EX) 2£ f + l ii bc k 



ds ' 2£i + l (2tt) 3 ^ 2 



n 



Y,Mk)Ky)(vly\\M(EX)\\v 9 .. 



(19) 



In order to test the accuracy of the THO basis to describe the continuum, we will calculate also the following 
magnitudes: 



The non-energy weighted sum rule: 



NEWSR : 



ds 



dB(EX) 
ds 



2£ ] + l ( W { )g - s - 



(20) 



with (r 2A ) g . s . = (<Pg. s .|r 2A |^ g . s .). 
• The energy-weighted sum rule: 



EWSR = 



dB(EX) 

dS ; (S — £g.s.) 



ds 



A(2A + 1) § L TT ( ^ )2(r2A " > --- 



(21) 
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• The energy- inverse weighted integral (or polarizability): 




1 dB(EX) 
(e - £ g . s .) de 



(22) 



Due to their respective weight factors, the EWSR and the polarizability are useful quantities to test the accuracy of 
the basis to describe high-energy and low-energy part of the spectrum, respectively. Note that the closed expression 
for the EWSR is only valid for angular momentum independent Hamiltonians. We note also that there is no closed 
expression for the polarizability, but in order to calculate this quantity with the desired accuracy, as well as the 
EWSR for angular momentum dependent Hamiltonians, one can directly evaluate (|22p using the exact continuum 
states integrated up to a sufficiently high excitation energy. 



As an illustration of the expressions derived in the preceding sections, we first consider the case of the p-n system 
with a central potential. Following [2lj | , the interaction between the proton and the neutron is parametrized in terms 
of the Poeschl-Teller potential, 



with V = 102.706 MeV and a = 0.9407 fin -1 . With these values, the ground state energy is 2.2245 MeV, in agreement 
with the experimental value. 

The oscillator length was chosen in order to minimize the ground state energy obtained upon diagonalization of 
the Hamiltonian in a small THO basis. This yields the value b = 1.5 fm. Once the value of b is fixed, the ratio -f/b 
determines the extension of the PS eigenstates; increasing the value of 7 reduces the radial extension and pushes the 
eigenvalues to higher energies. This is better seen in terms of the density of states, defined according to Eq. (j4j. This 
magnitude is plotted in Fig. Q] for the I = continuum, using a basis of N = 30 states, and three different choices 
of 7, namely, 7=1 fm 1 / 2 , 2.48 fm 1 / 2 and 5 fm 1 / 2 . It is seen that small values of 7 (which correspond to an extended 
THO basis in configuration space) produces a fine description of the continuum at low energy. This is useful, for 
example, to study Coulomb breakup. Increasing the value of 7 will decrease the density of states at low energies, 
that is compensated by an increase of the density at higher excitation energies. The most suitable choice for this 
parameter will depend on the problem at hand, depending on the energy region of interest. We emphasize, however, 
that the dependence on 7 is not critical and, in the applications shown here, different values of 7 converge to the same 
results for sufficiently large bases. For comparison, we include also in Fig. Q] the density obtained with a HO basis 
with TV = 30 states and b = 2.0 fm, which minimizes the ground-state energy for the HO basis. 

We next consider the scattering phase-shifts. These are displayed in Fig. [2] as a function of the relative p-n energy. 
The top, middle and bottom panels correspond to the s, p and d-waves. The solid line is the calculation using the 
asymptotics of the exact scattering states, whereas the circles represent the calculation obtained with Eq. ((?]), using 
a THO basis with N = 30 states. The LST was generated with the parameters b = 1.5 fm and 7 = 2.48 fm 1 / 2 . In 
the three cases, we find an excellent agreement between the exact and approximate phase-shifts in the whole energy 
range. 

Note that the calculated phase-shifts are consistent with the Levinson theorem (see e.g. [22)]), which establishes 
that the phase-shift at zero energy is given by Se(0) = mr, where n is the number of bound states for the partial wave 
£. So, since the s-wave supports a bound state (the deuteron), we have <5o(0) = n, whereas for £ — 1 and t — 2 we 
have <5i(0) = 5 2 (0) = 0. 

We now consider the electric transition probabilities, B(E1) and B(E2). These are shown in Fig. [3J The solid 
line corresponds to the calculation using the scattering states [Eq. fT3")) ]. the filled circles correspond to the discrete 
approximation using the THO basis [Eq. (IT71) ]. and the dashed line is the calculation obtained folding the discrete 
distribution with the continuum states [Eq. (1191) ]. Both the discrete and folded approximations show an excellent 
agreement with the exact distribution. We include also the calculation using the HO basis with N = 30 states and 
b = 2.0 fm (open circles and dot-dashed line). Both the B(E1) and B(E2) distributions depart significantly from 
the exact distributions. In addition, the HO basis produces a small density of states at low energy, which might be a 
drawback for scattering calculations. 

In Tables HI and HT1 we present the convergence of the ground state energy and the El and E2 sum rules with respect 
to the basis size. The last row lists the exact values obtained with the closed expressions of Eqs. (|2QI21j) . It is seen 



V. 



APPLICATION TO NUCLEAR STRUCTURE 



A. Application to the deuteron 
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FIG. 1. (Color online) Density of states for the I = continuum obtained with a THO basis with N = 30 states, defined 
according to Eq. Q for different choices of the 7 parameter. The dot-dashed curve is the density obtained with a HO basis 
with 30 states and b — 2 fm. See text for details. 




FIG. 2. (Color online) Phase-shifts for the deuteron system as a function of the relative p-n energy. The upper, middle and 
bottom panels correspond to £=0, 1 and 2 continuum states, respectively. See text for details. 
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FIG. 3. (Color online) Electric transition probabilities for the d=p+n system. The upper and bottom panels correspond to 
A = 1 and A = 2 transitions, respectively. See text for details. 



TABLE I. Convergence of the ground state energy and the total B(E1) and B(E2) transition probabilities for the deuteron 
case. 



N e gs (MeV) Total B(E1) (e 2 fm 2 ) Total B(E2) (e 2 fm 4 ) 





HO 


THO 


HO 


THO 


HO 


THO 


10 


-2.1570 


-2.2150 


0.81380 


0.85286 


14.6993 


17.9729 


20 


-2.2201 


-2.2245 


0.86200 


0.87129 


19.1257 


20.7938 


30 


-2.2237 


-2.2245 


0.86926 


0.87136 


20.2919 


20.8297 


40 


-2.2241 


-2.2245 


0.87079 


0.87136 


20.6422 


20.8297 


50 


-2.2242 


-2.2245 


0.87119 


0.87136 


20.7591 


20.8297 


60 


-2.2243 


-2.2245 


0.87130 


0.87136 


20.8015 


20.8297 


70 


-2.2243 


-2.2245 


0.87134 


0.87136 


20.8178 


20.8297 


Exact 




-2.2245 




0.87136 




20.8297 



that with a moderately small basis one obtains a very good convergence to the exact values. For comparison, in this 
Table we include also the calculations using the HO basis. From the quoted numbers, it is clear that the convergence 
rate is much faster for the THO basis. 

It is worth noting that, despite the simple Hamiltonian adopted in this work for the p — n system, the calculated 
polarizability is fully consistent with the experimental value a C xp = 0.61 ± 0.04, quoted in [23| . 
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TABLE II. Convergence of the polarizability and EWSR for the deuteron case. 



N 


HO 


a 
(fm 3 ) 


THO 


Ener; 

HO 


jy weighted B(E1) 
(e 2 fm 2 MeV) 

THO 


10 


0.413988 




0.515397 


8.605586 


8.102094 


20 


0.557520 




0.619350 


7.742321 


7.442755 


30 


0.599683 




0.620899 


7.526588 


7.437322 


40 


0.613204 




0.620922 


7.466005 


7.434683 


50 


0.617923 




0.620931 


7.446418 


7.433209 


60 


0.620109 




0.620936 


7.439133 


7.432331 


70 


0.620395 




0.620940 


7.435989 


7.431757 


Exact 




0.620953 






7.429937 



B. Application to 6 He 

We now consider a situation in which more complicated continuum structures are present, such as resonances. For 
this purpose, we take the 6 He nucleus, treated as a two-body system a+2n. Following (24[, the interaction between 
the two clusters is described with a Woods-Saxon shape, with R = 1.9 fm and a = 0.39 fm. For £ = states, the 
depth of this potential is adjusted to give the effective separation energy of 1.6 MeV between the two clusters. It was 
shown in Ref. [§4| that using this effective binding energy, instead of the two-neutron separation energy (SWi = 0.97 
MeV), provides a more realistic description of the ground state wave function. For I = 2, the inter-cluster potential 
is adjusted to yield a resonance at an excitation energy of E x — 1.8 MeV. For I = 1, we simply took the depth found 
for £ = 0. A THO basis with N = 50 states was used, and the LST was generated with the parameters b = 1 fm and 
7 = 1.89 fm 1 / 2 . 

The s-, p- and <i-wave phase shifts are displayed in Fig. @]as a function of the a+2n relative energy, e. For £ = 2, 
the energy scale has been restricted to the energy interval e = — 6 MeV in order to emphasize the region of the 
resonance. Again, we find a perfect agreement between the exact (solid line) and approximate (circles) phase-shifts in 
the whole energy range. Note that, in this case, the s-wave potential supports two bound states, the Pauli forbidden 
IS state and the 2S ground state. Consistently, the phase-shift at zero energy is given by S o (0) = 2ir. Analogously, for 
I = 1 the phase-shift tends to it, due to the presence of a (Pauli forbidden) bound state in this partial wave. Finally, 
for £ — 2 no bound states are supported by this potential and, therefore, ^(0) = 0. For the d-wave, the phase-shift 
crosses abruptly ir/2 at e — 0.20 MeV, reflecting the presence of a narrow 2 + resonance. Interestingly, this behavior 
is also observed in the THO basis, where there is a PS that appears exactly at the nominal energy of the resonance. 
It is then tempting to conclude that this PS will carry most of the character of the resonance, and in fact this is 
confirmed in Fig. [5j where we show the radial part of the exact scattering wave function, calculated at the energy of 
the resonance (solid line) along with the radial part of the PS eigenstate that appears at the energy of the resonance 
(dashed line). The former has been arbitrary normalized in order the two wave functions coincide at the maximum. 
It is seen that both wave functions are very similar up to very large distances. For comparison, we have included also 
the PS eigenstates associated to the eigenvalues just below (dotted-dashed line) and above (dotted line) the resonant 
one. They are very different from the scattering wave function at the resonance. In particular, it can be seen that a 
significant part of the norm of the resonant wave function is concentrated in the interior, as expected for a resonance, 
whereas for the non-resonant PS eigenstates, the probability in the interior is very small. This leads to the conclusion 
that in this case the character of the resonance is very well described by a single PS eigenstate. Indeed, if the basis is 
increased, the resonant character will be distributed among several PS. These results clearly show that a distinctive 
feature of the continuum such as the resonant structures are very well accounted for by the PS basis, despite its wrong 
asymptotic behaviour. 

In Fig. [5] we compare the El and E2 transition probabilities obtained from the scattering states by means of Eq. ([13")) 
(solid line), with the approximate distributions calculated with the THO basis. The circles correspond to the discrete 
expression of Eq. (|17p . whereas the dashed line is the smooth distribution obtained with Eq. ([IT?]) . Both the discrete 
and smooth distributions are in excellent agreement with the exact distribution. 

In Table Hill we display the convergence of the ground state energy and the El and E2 sum rules with respect to 
the basis size. As in the deuteron case, both observables converge very fast to their exact value, given by the sum 
rule of Eq. P0]). 

Finally, in Table ITVl we present the convergence of the integrated energy weighted B(E1) and the polarizability. For 
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FIG. 4. (Color online) Phase-shifts for the 6 He system as a function of the relative a+2n energy. The upper, middle and 
bottom panels correspond to £=0, 1 and 2 continuum states, respectively. See text for details. 



TABLE III. Convergence of the 


ground state enerj 


5y and the total B(E1) and B(E2) transition probabilities 


for the 6 He=a+2n 


system. 








N 


£gs 


Total B(EV) 


Total B{E2) 




(MeV) 


(e 2 fm 2 ) 


(e 2 fm 4 ) 


10 


-1.5913 


1.3538 


8.5524 


20 


-1.5999 


1.3854 


9.7471 


30 


-1.5999 


1.3855 


9.7471 


40 


-1.5999 


1.3855 


9.7471 


50 


-1.6000 


1.3855 


9.7471 


80 


-1.6000 


1.3855 


9.7471 


Exact 


-1.6000 


1.3855 


9.7471 



the former, we can not use the closed expression of Eq. (l21"j) , because the 6 He Hamiltonian depends on the angular 
momentum. The exact value listed in this Table corresponds to the explicit calculation using the scattering states 
up to a high excitation energy. The same holds for the polarizability, since for this observable there is no closed 
expression. 
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FIG. 6. (Color online) Electric transition probabilities for the 6 He=a+2n system. The upper and bottom panels correspond 
to A — 1 and X — 2 transitions, respectively. See text for details. 



VI. APPLICATIONS TO NUCLEAR REACTIONS 



A. Application to the reactions 6 He+ 12 C and 6 He+ 208 Pb at 240 MeV/nucleon 



The THO basis considered in this work is intended to provide a suitable discrete representation of the continuum 
spectrum of a loosely bound system, which can be useful for scattering calculations within the Continuum-Discretized 
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TABLE IV. Convergence of the polarizability and energy weighted B(El) for the 6 He=a+2n system. 
N a Total E. W. B(El) 





(fm 3 ) 


(e 2 fm 2 MeV) 


20 


1.8652 


6 5623 


30 


1.8746 


6.5464 


40 


1.8750 


6.5438 


50 


1.8752 


6.5425 


60 


1.8753 


6.5418 


70 


1.8753 


6.5412 


120 


1.8755 


6.5401 


Exact 


1.8756 


6.5393 



Coupled- Channels method Q. As a test case, we will apply the THO basis to the reactions 6 He+ 12 C and 6 He+ 208 Pb 
at 240 MeV/nucleon. These reactions were measured by Aumann et al. [25| at the GSI facility with the aim of 
extracting information on the 6 He nucleus. The breakup of 6 He on 208 Pb was already analysed using the CDCC 
method with the THO basis in Ref. |6j, showing an excellent agreement with the binning discretization method for 
the modulus of the breakup S-matrix. In this work, we extend the analysis of [6j in order to compare with the data 
of Ref. [25[ . In particular, we consider the exclusive breakup differential cross section as a function of the excitation 
energy of the projectile da/dE x . In Ref. (25|, this observable was obtained by reconstructing the kinematics of the 
6 He cm. from the measured momenta of the outgoing fragments ( 4 He+n+n) and integrating up to a laboratory 
scattering angle of 80 mrad. To obtain this observable in our calculations, we first construct the double differential 
cross section dcr/dQdE x from the breakup S-matrices. In principle, the breakup S'-matrix is a continuous function 
of the asymptotic momentum k. However, within a PS representation of the continuum, only discrete values of the 
S'-matrix are obtained, corresponding to the eigenvalues e n . A continuous breakup S-matrix can be obtained from 
the solution of the coupled equations following the procedure used in 0, [|| [26| , in which the discrete S-matrices are 
folded with the exact scattering states, similarly to what was done with the B(EX) distribution in Eq. (fl"9)) . i.e. 

S f:i (k) « f>*,(*OI^^(*»)> (24) 

71=1 

where S„^(fc„) are the discrete S— matrix elements resulting from the solution of the coupled-channels equations 
using a PS basis with N states. The subscripts i, n and / denote the channels {ip g . s .; Li, £i, J}, {ip^ ;Lf,£f, J} and 
{<pi] (k)', Lf, if, J}, respectively, where Li (Lf) is the initial (final) orbital angular momentum for the projectile-target 
relative motion, and J the total angular momentum of the system. 

We first consider the Pb target. For the 2n-a potential and the parameters of the THO basis we keep the same 
values of the preceding section. These parameters were found to provide an appropriate distribution of eigenstates for 
the 6 He+ 208 Pb reaction at 240 MeV/nucleon [(J. The a+ 208 Pb and 2n+ 208 Pb interactions are needed to generate the 
diagonal and non-diagonal coupling potentials of the CDCC equations. For the a+ 208 Pb interaction we have adopted 
the first optical potential used in Ref. [27j ■ For the 2n+ 208 Pb interaction, we used the following single- folding model: 

U(TL) = y , p„„(r n „){C/ n (R+^) + [/„(R-^)}dr n „ (25) 

where U n is the neutron- 208 Pb optical potential taken from the parametrization of Madland [28j], evaluated at the 
apropriate energy per nucleon, and p nn {r) is the neutron- neutron density distribution. The latter was calculated 
integrating the square of the three-body wave function of the 6 He nucleus along the 2n-a coordinate, i.e. 

p(r nn ) =r& n - a J |* 3fc (r„ n ,r 2 „_ Q )| 2 dr2n-a^n» , (26) 

where (r nn , T^n-a) is the three-body wave function and f2 n „ denotes the angular variables (0 nn ,4> nn ). In the 
present calculations, the function \l/ 3fc was taken from Ref. [24| . 

Very good convergence of the CDCC calculations was achieved with a basis of N = 30 states. In addition, we found 
that continuum states above 50 MeV have a negligible effect on the scattering observables, and hence these eigenstates 
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FIG. 7. (Color online) Angle-integrated breakup differential cross section, as a function of the excitation energy, for the reaction 
6 He+ 208 Pb at 240 MeV/nucleon. The dotted, dashed and dotted-dashed lines are the contribution of the s, p and d states. The 
thin solid line is the sum of these contributions and the thick solid line is the full calculation convoluted with the experimental 
resolution. The filled circles are the data from Ref. [25l ]. 



were removed from the coupled-channels calculation. This reduces the number of PS included in the CDCC equations 
to n s = n p = 14, n c i = 15 for s, p and d waves. The coupled equations were integrated up to 100 fm, and for a total 
angular momentum up to J max = 2000. 

In Fig. [7| we show the calculated energy differential breakup cross section, along with the GSI data. The dotted, 
dashed, and dotted-dashed lines are the separate contributions of the s, p and d continuum states. For this angular 
range (6\ a b < 80 mrad) the cross section is largely dominated by the coupling to the j = 1~ states due, mainly, to the 
strong dipole Coulomb interaction. 

The thin solid line in Fig. [7]is the sum of the s, p and d contributions. For a meaningful comparison with the data, 
this curve has to be convoluted with the experimental energy resolution, which we took from the same work (25j . The 
result of this folding is shown by the thick solid line. At low excitation energies, this calculation reproduces very well 
the shape and magnitude of the data. For excitation energies E x > 4 MeV, the calculation underpredicts the data. 
This discrepancy was also found in the semiclassical calculations reported in [25j . Note that the narrow peak due to 
the 2 + resonance in 6 He disappears in the folded calculation. 

The dominance of the dipole Coulomb couplings at these very small angles was used in Ref. (2f| to extract the 
dB(El) j 'de distribution, by comparing the measured differential cross section, da/dE x , with semiclassical calculations. 
Although our calculations confirm the dominance of the El couplings, we have found that nuclear potentials have 
a small but not negligible effect on this observable. In addition, starting from the same structure model for the 
6 He nucleus, our CDCC calculations show some departure from the semiclassical calculations, suggesting that the 
connection between the energy differential cross section and the underlying El probability is more complicated than 
suggested by the semiclassical approach. These results are potentially very interesting because they may affect the 
extracted dB(El) /de distribution from the cross section data. This analysis is beyond the scope of the present work 
and then we leave it for a separate publication. 

We consider now the 12 C target. The 2n-a interaction and the parameters of the LST where the same used in 
the Pb case. For the 6 He continuum states, £ — — 3 waves were included. In this case, we used a THO basis with 
N = 40 states. The ct+ 12 C potential was adopted from Ref. (29[. For the 2n+ 12 C interaction we used the same single- 
folding model as for the 208 Pb taking also the neutron-target potential, n-|- 12 C in this case, from the parametrization 
of Madland [28| . The coupled equations were integrated up to 100 fm, and for a total angular momentum up to 
Jmax = 200. In Fig. [5] we compare the CDCC calculations with the experimental data from (2f|. The meaning of 
the curves is the same as in Fig. [7] In this case, the low energy cross section is dominated by the population of the 
2 + continuum, with s and p waves contributing to the background. The contribution of the /-wave was found to be 
very small and hence it has not been included in this figure. As in the Pb case, the <i-wave cross section shows a 
narrow peak corresponding to the well known resonance at E x = 1.8 MeV. The width of the peak in our calculation 
is significantly smaller than the experimental width of this resonance (r exp s» 100 keV). This is a consequence of 
our simple two-body model adopted for the 6 He nucleus. Nevertheless, when the calculation is convoluted with the 
experimental resolution (thick solid line) it becomes very close to the data. Despite the simplicity of the structure 
model adopted in these calculations, these results show that the THO discretization method constitutes a useful 
method to describe accurately detailed structures of the continuum that may show up in scattering observables. 
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FIG. 8. (Color online) Energy breakup differential cross section, as a function of the excitation energy, for reaction 6 He+ 208 Pb 
at 240 MeV/nucleon. The dashed line is the CDCC calculation with the THO discretization method. The solid line is the 
result of convoluting the CDCC calculation with the experimental resolution. The filled circles are the data from Ref. [25[ . 



As a final remark on these results, we note t hat, at these energies, relativistic effects might affect the calculated 
observables. Some recent exploratory studies [13, HJ have shown that these effects produce an increase of about 10-15 
% on the calculated breakup cross section for the 8 B+ 208 Pb and n Be-|- 208 Pb reactions at 250 MeV/nucleon. These 
effects affect only the very small angles and are mostly due to the modification of the Coulomb potential. Therefore, 
we might expect a similar effect in our 6 He+ 208 Pb case. For the 6 He+ 12 C case, we do not expect these effects to be 
important, because dynamical relativistic corrections to the nuclear interaction were found by the same authors to 
be negligible at these energies. These corrections refer to dynamical effects exclusively. Relativistic kinematics effects 
were included in the referred works, as well as in our calculations, by using the appropriate relativistic momentum. 
In any case, the aim of our work is to show the ability of the THO basis to describe the continuum of a two-body 
system, and so the emphasis of our study is more on the description of the structure, rather than on the reaction 
mechanism. 



B. Application to the n Be+ 12 C reaction at 67 MeV/nucleon 

As a final example, we consider the scattering of the halo nucleus n Be on a carbon target. This reaction has been 
recently measured by Fukuda and collaborators at RIKEN [32[, in order to extract information on the 11 Be continuum 
by measuring neutron- 10 Be coincidences following the projectile breakup. The angle-integrated differential cross 
section as a function of the relative energy between the outgoing neutron and 10 Be, displays a structure dominated by 
a prominent resonance at E x = 1.78 MeV. This resonance was interpreted as a ^5/2 neutron coupled to the 10 Be ground 
state. A second bump was also observed at E x = 3.41 MeV, which was given a tentative assignment 3/2 + , with a small 
contribution of the d 3 / 2 wave coupled to 10 Be(0 + ) and a larger contribution of the 10 Be(2 + )(g)^2s 1 / 2 configuration. 
As in previous analyses of this reaction [33|, HH, we use a two-body model of the projectile, 11 Be= 10 Be(g.s.)+n, and 
hence those states based on the core excited states are absent from our model-space. 

In the present CDCC calculations, the neutron- 10 Be interaction was taken from (35[. This potential reproduces the 
separation energy of the ground state (l/2 + ) and first excited state (l/2~), and the position of the 5/2^ resonance, 
assuming the configurations 2s 1 / 2 , lpi/2 and d 5 / 2 , respectively. The 11 Be continuum was described with a THO basis 
with N = 25 states. The LST was generated with the parameters b = 2.4 fm and 7 = 3.6 fm 1 / 2 . Continuum states 
with configuration S\/ 2 , P1/2, P3/2, ^3/2 arL d d 5 / 2 were considered. After diagonalization of the projectile Hamiltonian 
in this THO basis, only those eigenstates with excitation energies below E < 20 MeV were included in the coupled- 
channels calculations, since the breakup cross section was found to be very small above this energy. This leaves about 
10-11 eigenstates for each partial wave. Following [34|, the n+ 12 C and 10 Be+ 12 C potentials were taken from [36| and 
(37j , respectively. The coupled-channels equations were integrated up to a matching radius of R = 90 fm and for total 
angular momenta up to J = 350. 

In Fig. [5] we compare the experimental [32| and calculated energy differential cross section as a function of the 
n- 10 Be relative energy. Both the data and the calculations correspond to the angular range 0° < # c .m. < 12°. The 
calculated contribution of each partial wave is shown. The symbols correspond to the contribution of specific pseudo- 
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FIG. 9. (Color online) Energy breakup differential cross section, as a function of the excitation energy, for reaction 11 Be+ 12 C 
at 67 MeV/nucleon. The separate contribution of each partial wave is shown. The dotted line is the full contribution, and the 
thick solid line is the folding of the latter with the instrumental resolution. The experimental data are from Ref. [3^ |. 



states, whereas the continuous lines are obtained convoluting the discrete ^-matrices with the exact continuum states. 
The low- lying continuum is dominated by the p 3 / 2 and d 5 / 2 waves, with the latter being responsible for the resonant 
peak at E x = 1.78 MeV. The dotted line is the sum of the different partial waves, and the thick solid line is the 
result of folding this full calculation with the instrumental resolution quoted in (32j . At energies close to the breakup 
threshold, the calculation overestimates the data. At energies above the resonance peak, the shape and magnitude of 
the data is well reproduced, with the exception of the broad peak due to the second resonance. We remind, however, 
that this resonance is believed to contain a significant contribution of 10 Bc(2 + ) and, therefore, it is not expected to 
be well described with our model space. Our results are very close to those obtained by Howell et al. [341 ]. including 
the overestimation of the data at low excitation energies. Since the calculations of that work used a continuum 
discretization in terms of energy bins, we attribute this discrepancy with the data to the choice of the interactions 
or to the restrictions of our three-body model, rather than to the method of discretization. As discussed in the 
6 He+ 12 C case, an advantage of the THO discretization over the standard binning method is the ability of describing 
fine structures of the continuum with a relatively small basis. For example, to describe d 5 / 2 resonance, the CDCC 
calculations of Ref. [HJ used 15 bins for e = 0.5 — 2 MeV, whereas in the present calculations about 10 PS are enough 
to describe the full energy region, including the narrow resonance. 

As in the 6 He case, for the d 5 / 2 continuum we get an eigenstate at e — 1.25 MeV, which is close to the nominal energy 
of the resonance (e = 1.27 MeV) and hence it is plausible to associate this eigenstate with the resonance structure. 
To corroborate this conclusion, in Fig. [TU] we compare the experimental angular distribution of the resonance region 
[32j with the angular distribution of the three eigenstates closer to the nominal resonance energy. As anticipated, the 
eigenstate at e = 1.25 MeV reproduces fairly well the shape and magnitude of the data, supporting our conclusion 
that this eigenstates carries most of the resonant character. It has to be borne in mind that, as the basis size is 
increased, the resonant character will be distributed among several eigenstates and hence this identification is not 
possible. In fact, for the present basis size (N — 25) there might be some mixing between the three eigenstates shown 
in Fig. [TUl 



VII. SUMMARY AND CONCLUSIONS 



In this work we have dealt with the problem of the study of the continuum properties of a weakly bound system 
in terms of basis of square-integrable functions, or pseudo-states (PS). The general idea of the PS method is to 
diagonalize the Hamiltonian of the two-body system in a truncated PS basis. The eigenstates with negative energies 
represent the bound states of the system, whereas those at positive energies are regarded as a finite and discrete 
representation of the continuum spectrum. Among the many possible choices of the PS basis, in this work we have 
made use of the Transformed Harmonic Oscillator (THO) basis proposed in [(| Hlj . in which the PS functions arc 
generated by applying a parametric Local Scale Transformation (LST) to the HO basis. The analytic form of the 
LST makes very simple the calculation of the PS basis. In addition, the radial extension of the basis and the energy 
distribution of the eigenvalues can be controlled through the parameters defining the LST. This permits to adapt the 
properties of the basis to the problem at hand. 
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FIG. 10. (Color online) Angular distribution for the E x = 1.78 MeV resonance. The circles are the data points from Ref. [32 
The curves are the angular distribution due to the three PS around the resonance. See text for details. 



In order to test the accuracy of the THO basis to represent the continuum, we have evaluated the scattering phase- 
shifts for the deuteron and 6 He systems, treated as two-body systems (p + n and a+2n) interacting with a simple 
central interaction. Since the THO states vanish at large distances, the phase-shifts have been evaluated using an 
integral formula, following the prescription of Hazi and Taylor [l2j . In both cases, we find an excellent agreement 
with the exact phase-shifts, obtained from the asymptotic part of the scattering states. Even the sharp resonance in 
the 6 He case is very well described with a small THO basis. As an additional test, we have evaluated the electric 
transition probabilities B(E\) and B(E2) for the same systems, finding again an excellent agreement with the results 
obtained with the scattering states. For this observable, a simple smoothing procedure has been proposed to provide 
a continuous distribution (dB(EX)/de) in terms of the discrete values obtained with the PS basis. 

Finally, we have presented CDCC calculations for the reactions 6 He+ 12 C and 6 He+ 208 Pb at 240 MeV/nucleon, 
and 11 Be+ 12 C at 67 MeV/nucleon for which experimental data exist [25|, [32[. For the n Be+ 12 C reaction, we have 
used a two-body model 10 Be(g.s.)+n of the projectile. In order to compare with the recent data of Fukuda et. [32] |. 
we have calculated the breakup differential cross section as a function of the neutron- 10 Be relative energy. Using a 
relatively small THO basis, we have been able to reproduce fairly well the data, including the narrow d 5 / 2 resonance 
at e — 1.27 MeV. Interestingly, one of the THO eigenstates appears at an energy very close to this energy and its 
associated differential angular cross section reproduces fairly well the experimental angular distribution obtained for 
the resonance region. 

For the 6 He reactions, we have used a simple two-body model (a+2n). Our calculations, which are parameter 
free, reproduce quantitatively and qualitatively the experimental energy differential cross sections reported in [25[ 
both for the heavy target ( 208 Pb) and for the light target ( 12 C). Furthermore, these calculations indicate that, for 
the 6 He+ 208 Pb reaction at excitation energies below ^4 MeV, break-up cross sections are Coulomb dominated, with 
monopole and quadrupole components contributing only about 6%. However, for 6 He+ 12 C, the dominant component 
is the quadrupole, so that for excitation energies below ~2.5 MeV, dipole and monopole components contribute about 
24%. 

The role of nuclear forces and higher order effects has been investigated by comparing our full coupled channels 
calculation with Coulomb calculations using the equivalent photon model. Differences as large as 28% have been 
found, indicating the need of performing continuum-discretized coupled-channels calculations to extract structure 
information from break-up reaction data, at least in this energy regime. 

The agreement between theory and experiment is very encouraging, given the simplicity of the dineutron model used 
in the present calculations. One has to bear in mind, however, that some of the details of the breakup distributions 
might be hidden due to the energy resolution of the experiment. New measurements with better energy resolution 
will be useful to test more stringently the break-up distributions at energies closer to the threshold. An accurate 
description of these reactions will require a realistic three-body model to describe the 6 He nucleus. In this respect, it 
is worth noting that the THO method used in this work can be generalized to three-body problems. This can provide 
a useful yet simple method to study continuum structures (e.g. resonances) in nuclei with a three-body structure 
( 9,14 Be, 6 He, 11 Li, 8 B, etc) as well as reactions involving these nuclei. A similar approach proposed very recently, 
making use of a different PS basis, has been found to provide very promising results (38l. [39j. 
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